%% dx = (0.125,0.125)
solution = [
1 1.125 1.25 1.375 1.5 1.625 1.75 1.875 2 
1.125 1.25 1.375 1.5 1.625 1.75 1.875 2 2.125 
1.25 1.375 1.5 0 0 0 2 2.125 2.25 
1.375 1.5 1.625 0 0 0 2.125 2.25 2.375 
1.5 1.625 1.75 0 0 0 2.25 2.375 2.5 
1.625 1.75 1.875 2 2.125 2.25 2.375 2.5 2.625 
1.75 1.875 2 2.125 2.25 2.375 2.5 2.625 2.75 
1.875 2 2.125 2.25 2.375 2.5 2.625 2.75 2.875 
2 2.125 2.25 2.375 2.5 2.625 2.75 2.875 3 
];
exact = [
1 1.125 1.25 1.375 1.5 1.625 1.75 1.875 2 
1.125 1.25 1.375 1.5 1.625 1.75 1.875 2 2.125 
1.25 1.375 1.5 0 0 0 2 2.125 2.25 
1.375 1.5 1.625 0 0 0 2.125 2.25 2.375 
1.5 1.625 1.75 0 0 0 2.25 2.375 2.5 
1.625 1.75 1.875 2 2.125 2.25 2.375 2.5 2.625 
1.75 1.875 2 2.125 2.25 2.375 2.5 2.625 2.75 
1.875 2 2.125 2.25 2.375 2.5 2.625 2.75 2.875 
2 2.125 2.25 2.375 2.5 2.625 2.75 2.875 3 
];
err = [
0 1.33227e-15 -3.55271e-15 3.33067e-15 1.37668e-14 -3.10862e-15 -3.17524e-14 1.46549e-14 0 
4.44089e-16 2.22045e-15 1.55431e-15 2.66454e-15 3.10862e-15 -6.43929e-15 -1.15463e-14 0 0 
5.77316e-15 5.32907e-15 4.88498e-15 0 0 0 -7.77156e-15 -2.66454e-15 0 
0 8.43769e-15 9.76996e-15 0 0 0 -2.66454e-15 -1.33227e-15 0 
4.44089e-16 1.84297e-14 1.9762e-14 0 0 0 3.55271e-15 1.33227e-15 0 
1.38334e-13 4.52971e-14 2.84217e-14 2.30926e-14 1.77636e-14 1.19904e-14 7.10543e-15 3.10862e-15 0 
-1.02141e-13 -3.9968e-15 2.57572e-14 2.35367e-14 1.82077e-14 1.77636e-14 9.76996e-15 3.9968e-15 0 
9.76996e-15 1.5099e-14 5.4623e-14 2.70894e-14 1.42109e-14 2.93099e-14 8.43769e-15 3.10862e-15 0 
0 0 1.5099e-13 1.64313e-14 -1.77636e-14 7.72715e-14 -8.43769e-15 0 0 
];
